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Abstract 

C-function MADS-box transcription factors belong to the AGAMOUS (/AG) lineage and specify both stamen and carpel 
Identity and floral meristem determinacy. In core eudicots, the AG lineage is further divided into two branches, the euAG 
and RLE lineages. Functional analyses across flowering plants strongly support the idea that duplicated AG lineage genes 
have different degrees of subfunctionalization of the C-function. The legume Medicago truncatula contains three C-lineage 
genes in its genome: two euAG genes {MtAGa and MtAGb) and one PLENA-Wke gene {MtSHP). This species is therefore a 
good experimental system to study the effects of gene duplication within the AG subfamily. We have studied the respective 
functions of each euAG genes in M. truncatula employing expression analyses and reverse genetic approaches. Our results 
show that the M. truncatula euAG- and PLENA-Wke genes are an example of subfunctionalization as a result of a change in 
expression pattern. IVltAGa and IVltAGb are the only genes showing a full C-function activity, concomitant with their 
ancestral expression profile, early in the floral meristem, and in the third and fourth floral whorls during floral development. 
In contrast, MtSHP expression appears late during floral development suggesting it does not contribute significantly to the 
C-function. Furthermore, the redundant MtAGa and MtAGb paralogs have been retained which provides the overall dosage 
required to specify the C-function in M. truncatula. 
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introduction 

Genetic regulation of flower development has been subject of 
study over the last decades, particularly in the model species 
Ambidopsis thaliana and Antirrhinum majm. These studies 
provided a general understanding of floral organ development in 
higher plants and led to the proposal of the ABCDE model, which 
postulates that floral organ identity in each whorl is defined by five 
functions named A, B, C, D and E (for review [1,2,3,4]). In 
particular, C-function is required to promote stamen and carpel 
identity, to establish the determinate nature of the floral meristem 
[5] and also to repress A-function in the third and fourth whorls 
[6]. Cloning of ABCDE organ identity genes in Ambidopsis 
showed that most of them encode MADS-box transcription 
factors. Studies of MADS-box genes in higher eudicotyledoneous 



flowering plants show that they are key regulators of flower 
development. 

In Arabidopsis and Antirrhinum, the C-function is essentially 
represented by a single gene respectively, AGAMOUS (AG) [5,7] 
and PLENA {PLE) [8]. Additional C-function genes have been 
identified: SHATTERPROOFl {SHPl) and SHATTER- 
PR00F2 (SHP2) [9,10] genes in Arabidopsis and FARINELLI 
(FAR) [1 1] in Antirrhinum. 

Phylogenetic studies using a large data set oi AG-Yikt sequences 
show that AG and PLE actually represent paralogous lineages 
derived from a duplication in a common ancestor early in the 
history of the core eudicots. This duplication gave rise to the eu^G 
lineage which includes AG and EAR, and the PLENA lineage 
[PLE), where SHP1/SHP2 and PLENA are placed 
[12,13,14,15]. An even more ancient duplication occurred before 
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the radiation of extant angiosperms, producing the C Kneage {AG 
lineage) and the ovule-specific D lineage {AGLll lineage) [12,15]. 

AG, FAR and PLE all display very similar expression patterns 
during the developing male and female reproductive organs 
[5,8,1 1]. Mutations of AG and PLENA in their respective species 
producx; identical phenotypes characterized by developing flowers 
with petals in whorl 3 instead of stamens, and sepals in whorl 4 
instead of carpels. In addition, the floral meristem is not 
determinate [7,8]. However, far mutation only afiects male 
reproductive organs, causing partial male sterility [11]. Mean- 
while, SHP genes are only expressed in the ovules, in the 
developing pistil and show almost identical expression patterns in 
developing Arabidopsis fruit [16,17,18]. They function redundant- 
ly in style and stigma development, notably in the fusion of the 
carpel [19], and in seedpod shattering [10]. SHP and FAR genes 
do not contribute significantly to the C-function although FAR is 
expressed at early stages of flower development. 

These studies showed a random evolutionary trajectory for gene 
functions after a duplication event. AG and PLE, difierent 
members of a duplicated gene pair, retained the primary homeotic 
functions in different lineages. Their respective orthologs have 
taken completely new roles in fruit dehiscence, as the case of 
SIIPl / SHP2 genes, or contribute redundantly in male repro- 
ductive development, as the case oiFAR. Representatives of both 
eu^G and PLE lineages have been identified in different core 
eudicot species and several functional analyses are available for the 
paralogs from petunia, tomato, and Nicotiana benthamiana 
[20,21,22,23]. 

In addition to the ancient duplication events observed in the 
AGAMOUS subfamUy, more recent duplication events appear to 
have occurred in the euAG lineage [12,15,24]. There have been 
reported two members of the eu^G lineage in some species, as for 
example Populus trichocarpa [25], Cucumis salivus [26], Gerhera 
hybrida [27], and in the legume species Lotus japonicus [28], 
Pisum sativum and Medicago truncatula (this study). However, in 
most cases, functional analyses are lacking. Both phylogenetic and 
functional data of gene lineages are very important to understand 
the evolution of gene famflies [12,29]. 

Medicago truncatula {Mt) has three C-hneage genes: two euAG 
genes (MtAGa and MtAGb) and one PLENA-Vike gene {MtSHP; 
[30]). To gain insight into the specific contribution of the eu^G 
and PZ,£-like paralogous genes in the control of the C-function, 
we compared the expression patterns of the three C-hneage genes 
during flower development. We have also characterized MtAGa 
and MlAGh loss-of-func tioii mutants and plants where both genes 
have been simultaneously sflenced. The particular capacity of 
these two genes to promote reproductive organ identity has been 
tested by ectopicaUy expressing them in Arabidopsis. Our results 
indicate that the members of evAG and PLE lineages in M. 
truncatula are subfunctionalized, where the C-function is only 
promoted by the eu^G paralogs. They largely overlap in function 
but the overaU dosage of both gene products is critical to promote 
complete stamen and carpel identity and floral meristem 
determinacy in M. truncatula. 

Material and Methods 

Plant material and growth conditions 

Medicago truncatula cv. Jemalong lines A17, SA1335 and 
R108, and Arabidopsis thaliana cv. Columbia plants were used in 
this study. Plants were grown in the greenhouse, at 22°C (day) and 
18°C (night) with a 16 h light/8 h dark photoperiod, in soiI:sand 
(3:1) irrigated with Hoagland N".l solution supplemented with 
oligoelements [31]. 



The mtaga mutant allele (previously matag-2) was isolated in a 
previous screening [32] and homozygous plants were used in this 
study. 

Identification of Tnti insertion sites in MtAGb and co- 
segregation test 

The M. truncatula population used for the screening of mutants 
was described in detail [33,34,35,36] (http:/ /bioinfo4.noble.org/ 
mutant/). The mtagb aUele was identified by PCR screening of a 
segregating population of approximately 10,000 independent Unes, 
using primers annealing to the MtAGb sequence (AGb-F, Table 
S2) in combination with primers annealing to the LTR borders of 
the Tnll retroelement (Tntl-F; Table S2 and Figure S3). We 
identified a line (NF4908) with an insertion of the retroelement 
located in the first intron at 277 bp of the starting codon 
(Figure lA). The Rl plants were genotyped by PCR using the 
Tntl-F primer in combination with the gene-specific primers 
MtAGb-F and MtAGb-Rgenomic (Table S2 and Figure S3B). 
Approximately 70% of the plants showed a mutant phenotype and 
co-segregated with the Tntl insertion. 

Isolation and sequence analysis 

MtAGa and MtAGb cDNAs were isolated from a library of M. 
truncatula All inflorc-sccnce apices [37], using the MADS-box 
fragment of the M. truncatula PISTILLATA gene as a probe. 
Sequence alignments and similarity comparisons of the inferred 
proteins were performed using Align and ClustalW tools [38] . The 
deduced amino acid sequences were aUgned using ClustalW and 
further refined by hand. Genomic sequences search was 
performed using BLAST [39] . 

The second intron sequences were obtained by PCR using 
genomic DNA and the primers MtAGa-intron-F and MtAGa- 
intron-R for MtAGa and MtSHP-intron- F and MtSHP-intron-R 
for MtSHP (Table S2). To predict the probabflity of the presence 
of LFY binding sites in the intron sequences, we used the 
Morpheus webpage facility (http://biodev.cea.fr/morpheus/ 
Default, aspx) with the LFY matrix. The score tool permits 
localization of the best transcription factor binding sites in DNA 
fragments. 

Phylogenetic tree 

The phylogenetic tree was inferred by Neighbor-Joining using 
Poisson-Corrected amino acid distances. Reliabflity of internal 
nodes was assessed using bootstrap witii 10000 pseudo-replicates. 
Tree inference was conducted using MEGA version 4 [40] . The 
data set comprised 28 previously reported C- and D-class genes 
obtained from GenBank and the two new sequences that we 
isolated {MtAGa and MtAGb). The tree was rooted using the 
Arabidopsis D-class gene SEEDSTICK {STK) sctiuciice. All 
sequences used in this analysis, with their respective species and 
accession numbers, are included in Table SI. 

Southern blot hybridization 

Plant genomic DNA was extracted from leaves as described by 
[41]. Ten micrograms of DNA were digested with restriction 
enzymes, separated on 0.7% Tris-borate EDTA IX agarose gels 
overnight at IV/ cm and transferred to a nylon membrane. 
Southern blot hybridization was performed by standard methods 
at 65°C. A 241 bp fragment oi MtAGa cDNA (positions 572-813) 
and a 215 bp fragment of the MtAGb cDNA (positions 558-773) 
were amplified with the specific primers MtAGadir, MtAGarev, 
MtAGbdir and MtAGbrevl (Table S2) and used as probes. 
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Figure 1. Two euAGAMOUS genes in Medicago truncatula. (A) Gene structure of MtAGa and MtAGb. Coding sequences are represented as boxes 
and introns as dotted lines. Black triangles localize the position of the Tntl insertions present in the mutant lines mtaga (NF13380) and mtagb 
(NF4908) used in this study. (B) Neighbor-Joining Tree of euAG and PLENA homologs from a selection of diverse species. The numbers next to the 
nodes refer to bootstrap values from 10000 pseudo-replicates. (C) Distribution of putative LFY binding sites in the first intron of MtAGa, MtAGb and 
MtSHP genes as identified with the use of a position-specific scoring matrix using a cutoff value of —20. 
doi:1 0.1 371 /journal.pone.01 03770.g001 



Expression analyses 

Total RNA was isolated from frozen plant material using the 
RNeasy Plant mini Kit (Qiagen, Germany) according to the 
manufacturer's instructions. 



For Northern blot analysis total RNA (15 |Ig) from frozen 
leaves, roots, stems, flowers and poUinated ovaries (young fruits) 
was used. RNA electrophoresis was carried out in formaldehyde- 
agarose gels, transferred to Hybond N^ membranes (Amersham 
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Biosciences, USA), and hybridized with P-labeled probes under 
standard conditions. The probes were generated from the same 
gene fragments used for Southern blot analysis. 

For Real Time RT-PCR analysis total RNA was treated with 
rDNasel of the DNa.se Treatment and Removal Kit (Ambion, Life 
Technologies, USA). For first-strand synthesis, total RNA (1 |j,g) 
was reverse-transcribed in a 20 (tl reaction mixture using the 
PrimerScript 1*' strand cDNA Synthesis Kit (Takara, Japan). One 
microliter of RT reaction was used for a Real Time RT-PCR 
analysis with 300 nM of each primer mixed with the Power SYBR 
Green PGR Master Mix (Applied Biosystems) according to the 
manufacturer's instructions. The reaction was carried out into 96 
well-optical reaction plates using an ABI PRISM 7500 Sequence 
Detection System and appropriate software (Applied Biosystems). 
The relative levels were determined by the 2 "^'^^ Method [42]. 
To normalize the variance among samples, Secret Agent [O-linked 
N-acetyl glucosamine transferase: TG77416; [43]) was used as an 
endogenous control. All reactions were performed by triplicate 
using a biological replicate for each sample. Primers were designed 
using PRIMER EXPRESS software (Applied Biosystems, USA) 
using default parameters and are listed in Table S2. 

RNA in situ hybridization was performed on 8 |iM paraSin 
sections of M. truncatula inflorescences as described by [16], using 
digoxigenin-labeUed probes. The RNA sense and antisense probes 
were generated from the same gene fragments used for Southern 
blot analysis. Both fragments were cloned into the pGEM T-easy 
vector (Promega, USA) and the probes were synthesized using SPG 
or T7 RNA polymerases. 

Virus Induced Gene Silencing in Medicago truncatula 

pCAPE 1 and pC APE2 derivatives were used as vectors for gene 
silencing [44]. Two DNA fragments from the 3' region of the 
MtAGa (310 bp) and MtAGh (338 bp) genes were obtained by 
PGR using primers (MtAGaVIGSdir, MtAGbVIGSdir, MtAGa- 
VIGS2rev, MtAGbVIGS2rev) that added restriction sites to both 
ends of the fragments (Table S2). The amplicons were cloned into 
pGEM T-easy (Promega, USA), digested using the appropriated 
restriction enzymes and cloned into a similarly digested pGAPE2 
vector [44]. The resulting plasmid (pGAPE2-Mi/4Ga6) was 
confirmed by sequencing before being introduced into Agrobac- 
terium lumefaciens strain C58/pMP9(). Agrohaclerium inoculation 
of M. truncatula leaves was performed as described by [44]. 

Generation of transgenic RNA! plants 

Transformation of M. truncatula R108 was performed as 
described previously [33]. The 35S::RNAi-M</lG construct was 
performed using a 215 bp fragment from MtAGb (positions 557- 
772 from the ATG codon), amplified using primers MtAGb- 
RNAiD and MtAGb-RNAiR (Table S2) that incorporate two 
restriction sites that are used for cloning into the pHANNIBAL 
vector [45]. 

Arabidopsis transformation 

MtAGa and MtAGh cDNA fragments were amplified using 
AGaSBXdir, AGaSBXrev, AGbSBXdir and AGbSBXrev primers 
(Table S2) and cloned into the pBINJITGO vector [46], a pBIN19 
derivative (Glontech, Palo Alto, GA, USA), which placed gene 
transcription under the control of a tandem repeat of the GaMV 
35S promoter. Arabidopsis plants were transformed by floral 
dipping according to standard procedures [47] and selected on 
kanamycin. 



Photography, microscopy and cryo-SEM 

Light photographs were made with a stereomicroscope Leica 
MZ16F attached to a DFG300 FX camera (Leica Microsystems, 
Germany). 

After in situ hybridizations, sections were observed and 
photographed with a Nikon Eclipse E-600 microscope equipped 
with a digital camera. 

For cr^'o-SEM, samples were frozen in slush nitrogen and 
attached to the specimen holder of a Cr\oTrans 1500 Cryo- 
Preparation System (Oxford Instruments, UK) interfaced with a 
JEOLJSM-5410 scanning electron microscope. The samples were 
then transferred from cryostage to the microscope sample stage, 
where the condensed surface water was sublimed by controlled 
warming to -85°C. Afterwards, the sample was transferred again 
to the cryostage in order to gold coat it by sputtering. Finally the 
sample was put back on the microscope sample stage to be viewed 
at an accelerating voltage of 15 KeV. 

Results 

Identification of two euAGAMOUS genes in Medicago 
truncatula 

To identify MADS-box genes involved in flower development, 
we screened a cDNA library of M. truncatula floral apices using a 
set of MADS-box fragments from different species as a probe 
[48,49]. Among the isolated clones, three corresponded to full- 
length sequences that present significant similarity to genes of the 
AGAMOUS subfamily. One of these clones presents high 
sequence similarity with the Arabidopsis SHPl and SHP2 genes 
[17,18,50] and with the Antirrhinum PLENA gene {PLE; [8]) and 
corresponds to the MtSHP gene [30]. The other two clones 
isolated present high sequence similarity with the G-lineage genes 
AGAMOUS ot Arabidopsis [5] and FARINELLI of Antirrhinum 
[1 1] and have been named MtAGa and MtAGh. The MtAGa 
clone is 1208 bp long, with an ORE of 780 bp and a deduced 
protein of 260 amino acids. The MtAGh clone is 1099 bp long, 
with an ORF of 732 bp and the deduced protein has 244 amino 
acids. MtAGa shows 81 % amino acid identity with MtAGb, being 
more similar in the N-terminal domains (98% in the MADS 
domain, 76% in the I region and 80% in the K domain) than in 
the G-terminal region (68% of identity). MtAGa and MtAGb show 
respectively 71 and 68% amino acid identity with the euAG clade 
protein FAR and either one 67% amino acid identity with AG 
(Figure SI). 

The MtAGb genomic sequence was found in two BAG clones: 
mth2-30e7 (GenBank AC 137837.4) and mth2-76i7 (GenBank 
AC 153460.24). This gene is organized in seven exons and six 
introns (GenBank KJ470634; Figure lA). The genomic sequence 
of the MtAGa gene was not available in the databases and was 
obtained in this study (GenBank KJ47()633). MtAGa gene is also 
organized in seven exons and six introns (Figure lA). Both genes 
are present as single copy in the genome as confirmed by Southern 
blot analysis (Figure S2A). Genomic sequences from MtAGa and 
MtAGb genes were used as BLASTN queries, and displayed using 
the Chromosome Visualization Tool (CViT, http://www. 
medicagohapmap.org/tools/blastform). We found that MtAGh 
gene is located on chromosome 8 (Figure S2B). However, no 
location was obtained for the MtAGa gene. 

Phylogenetic analyses showed that MtAGa and MtAGb are 
relatively recent paralogs within the euAG subclade (Figure IB). 
Other model l(■gum(^s such as Pisum, Lotus and Glycine used in 
this analyses, also displayed two paralogs in this subclade. In most 
of the studied plants, the euAG subclade is represented by one 
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single gene, such usAG for Arabidopsis, FAR for Anthirrinum or 
TAGl for tomato. 

Genes from the eu^G- and -PLi?-subclades present a well 
conserved gene structure and characteristic large introns essential 
for their correct expression pattern [51,52]. InM. truncatula, both 
MtAG genes show a large first intron (4607 bp in MtAGa and 
3007 bp in MtAGb). We searched in these sequences for the 
presence of putative LFY binding sites that are characteristic of C- 
function genes. We included in this study the MtSHP gene and 
then we have sequenced its first intron of 4281 bp (GenBank 
KJ470635). We used a bioinformatics tool available at the 
Morpheus webpage facility (http://biodev.cea.fr/morpheus/ 
Default.aspx) which examines the distribution of individual LFY 
binding sites. The presence or absence of LFY binding sites 
identified by this approach appears to be helpful in predicting to 
what extent a C-clade gene acts as a true C-function gene [53]. 
Our results show different binding site landscapes for the three 
genes analysed (Figure IC). In the case of the MtAGb first intron, 
there is a main single binding site of very high score similar to the 
profile found for the PLE gene [53]. This site corresponds to two 
fully intact CCAAT-boxes. In the first intron of MtAGa, lower 
affinity sites are present but the binding seems to be compensated 
through tlie ai tion of several nearby sites as reported for the sites 
present in the- scxond intron of^G (Figure IC). The prediction of 
binding sites for LFY in the first intron oi MlSHP gave only low 
score values (Figure IC) indicating a low probability to be 
regulated by LFY at early stages of floral development. 

Globally, these data suggest that the duplicated e^^AG genes 
MtAGa and MtAGb might contribute to C-function specification 
in M. truncatula and play a role during early stages of flower 
development. 



MtAGa 



L Fl Fr S R 




MtAGb 

L Fl Fr S R 




Expression patterns comparison of three M. truncatula 
C-lineage genes during floral development 

The expression patterns of MtAGa, MtAGb and MtSHP were 
analysed by Northern blot in different plant tissues and the three 
genes are exclusively expressed in floral and young fruit tissues 
(Figure 2). We performed detailed in situ hybridization experi- 
ments to show the distribution of MlAGa, MlAGb and MtSHP 
mRNAs during floral development (Figure 3). Both MtAGa and 
MtAGb transcripts began to accumulate at stage 2 of flower 
development. At this stage, MtAGb accumulates in the centre of 
the floral primordia while MtAGa signal was detected throughout 
the floral meristem (Figure 3F and 3A). At stage 4, MtAGb 
transcript was located in the region of the common primordia that 
will give rise to the stamens and also in the central part of the floral 
apex where the carpel is developing (Figure 3G). However, 
MtAGa expression is still observed on the whole floral meristem, 
including petal and sepal primordia (Figure 3B). From stage 5, 
expression of both paralogs was distributed uniformly in whorls 3 
and 4 (Figure 3C-E and 3H-J), although hybridization signal for 
MtAGb was stronger on the abaxial region of the carpel. In later 
stages, expression of both transcripts was observed in the 
developing ovules, in the distal region of the carpel and on the 
filament of the anthers (Figure 3E and 3J). In contrast, MtSHP 
mRNA began to accumulate late on flower development and can 
be detected at stage 6 in the inner cells of the developing carpel 
(Figure 3M). Since late stage 7, MtSHP expression is exclusively 
detected in the ovules (Figure 3N). 

We conclude that the expression of the two M. truncatula eu^G 
genes remained largely restricted to the developing male and 
female organs from early stages and throughout flower develop- 
ment. MtAGa and MtAGb showed nearly identical spatial 
expression patterns from stage 4 to late developmental stages. 



MtSHP 

L Fl Fr S R 




Figure 2. Expression patterns of MtAGa, MtAGb and MtSHP 
genes in various plant tissues of M. truncatula. Northern blot 
analyses were performed using total RNA prepared from leaves (L), 
flowers (Fl), young fruits (Fr), stems (S) and roots (R). 
doi:1 0.1 371/journal.pone.01 03770.g002 

Differences were mainly observed at early stages, where MtAGa 
expression appeared to be wider spread in the floral meristem than 
MlAGb. Essentially, the duplicated MtAGa and MtAGb genes 
showed similar spatial and temporal expression patterns to those 
described for C-function genes from the eu^G subclade such as 
AG and PLE. In contrast, MtSHP was only expressed in the 
developing carp(-l and in the ovules, similar to the pattern 
described for Arabidopsis SHP genes that do not significantly 
contribute to C-function. 
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Figure 3. Expression pattern of MtAGa, MtAGb and MtSHP genes during early floral development. In situ localization of MtAGa (A-E), 
MtAGb (F-J) and MtSHP (K-N) transcripts in M. truncatula wild-type flower buds. Developmental stages were defined according to [74]. MtAGa 
transcripts are localized in the whole floral meristem at stages 2 (A) and 4 (B). At stage 5 (C) MtAGa mRNA is detected in stamen and carpel primordia. 
At stage 7 (D-E) expression locates in stamens, carpel and the developing ovules. MtAGb transcripts are localized in the center of the floral meristem 
at stage 2 (F). At stage 4 (G) expression locates in the carpel primordia and the half of the common primordia that will give rise to the stamens. At 
stage 5 (H) MtAGb mRNA is detected in stamen and carpel primordia. At stage 7 (l-J) MtAGb mRNA is detected in stamens, carpel and developing 
ovules. MtSHP transcripts are not detected at stages 2 (K) and 5 (L). At stage 6 MtSHP transcripts are localized in the inner cells of the developing 
carpel (M). At stage 7 (N) expression can be detected in the developing ovules. F: floral meristem; S: sepal; C: carpel; CP: common primordia; St: 
stamen; P: petal; Fi: filament; A: anther; Ov: ovule. 
doi:1 0.1 371 /journal.pone.01 03770.g003 



MtAGa and MtAGb loss-of-function analyses 

To investigate the specific contribution oi MlAGa and MtAGb 
genes in the M. truncatula floral development, we looked for 
retrotransposon insertion mutants [32]. These mutants were 
available for both genes, each containing a single Tntl element 
inserted in the N-terminal part of the gene (Figure lA). mtagh 
Tntl insertion was identified as described in the experimental 
procedures and it was located in the frrst intron at 277 bp from the 
start codon (Figure 1 A), mtaga (previously mtag-Z), was isolated as 
a case study to demonstrate the utility of the reverse genetic 
platform in the model legume Medicago truncatula [32]. Using 
quantitative RT-PCR we have confirmed that MtAGa or MtAGh 
transcripts were nearly undetectable in the corresponding homo- 
zygous mutant plants without affecting the expression levels of the 
corresponding MtAG paralog (Figure S4). tntaga and mtagh 
mutants exhibit very similar floral phenotypes: flowers are wild- 
type in appearance, showing only mild developmental defects on 
die third and fourth-whorl organs (Figure 4). Incomplete fusion of 
the staminal tube was observed and occasionally, some petaloid 
prolongation in the anther tip (Figure 4B and 4C). Whorl 4 is 
frequently composed by multiple carpels (2-3) or by modified 
carpels with stigmatic protuberances and exposed ovules (Fig- 
ure 4B and 4C). Flowers were generally sterile, although 
occasionally a few seeds were produced. 

To investigate the effect of the combined loss-of-function of 
MtAGa and MtAGb, we performed virus induced gene silencing 



(VIGS) experiments. We used a VIGS vector based on PEBV {Pea 
early browning virus) initially developed to induce gene silencing 
in Pisum [44] . The silencing efficiency of this vector is lower in M. 
truncatula cultivars than in P. sativum [54]. We generated the 
pCAPE2-Mt4 Ga6 construct to simultaneously silence both genes 
(see Materials and Methods). Twenty one plants were inoculated 
with the -pCAPYJi-MtAGah construct and we analysed the 
phenotype of 150 flowers. 10% of the MtAGab-VIGS flowers 
showed homeotic transformations in whorls 3 or 4. Only 1% of the 
MtAGab-VIGS flowers showed a near complete loss-of-C-function 
phenotype (Figure 4 D). In whorl 3 some anthers are completely 
converted into petal-like organs, and in the fourth whorl, sepal-like 
structures arise instead of carpels (Figure 4 D). These flowers also 
showed indeterminacy of the floral meristem revealed by the 
presence of multiple petaloid or sepaloid concentric structures in 
the centre (Figure 4D, arrows). This phenotype suggests that both 
MtAG genes were strongly down-regulated but, because of floral 
tissue limitation, we could not determine the expression levels of 
the targeted genes. 

In parallel, we also developed transgenic M. truncatula lines in 
which MtAG genes were down regulated by RNA interference 
(RNAi) and three independent lines were obtained. Only one of 
these lines (line 5.7) showed floral phenotypes similar to those of 
single mtag mutants (Figure S5). These flowers showed clear 
morphological aberrations consisting of mild homeotic transfor- 
mation of stamens to petals and multiple carpels with exposed 
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Figure 4. Floral phenotypes of M. truncatula mtagb and mtaga mutants and MtAGab-VIGS silenced plants. (A) Dissected wild-type M. 
truncatula flower showing the four floral whorls (Wl to W4). Floral phenotypes are similar in mtaga. (B) and mtagb (C) mutants. Mutant flowers (left) 
were opened to show the inner whorls. In W3 the staminal tubes are unfused and occasionally petaloid structures appear (arrow) replacing the 
anthers. Carpels in W4 present stigmatic protuberances (arrow) or multiple unfused carpels and exposed ovules (boxed). (D) /MMGob-VIGS flower with 
severe homeotic transformation of stamens into petals and carpels into sepaloid structures (arrows). Bars indicate 1 mm. 
doi:1 0.1 371/journal.pone.01 03770.g004 



ovules. We analyzed the expression of both MtAGa and MtAGh 
genes in the three RNAi lines. These lines showed a difTerent 
degree of silencing for each targeted gene (Figure S5J). In silenced 
line 1.3, we found MtAGb levels considerably reduced (~I5% of 
wild-type level), whUe MtAGa levels are somewhat decreased 
(~9()"/o of wild-type level). Strong silencing of one paralogous gene 
did not produce homeotic floral defects. However, floral defects in 
line 5.7 correlate with simultaneous downregulation of both genes 
(-40% of wild-type level) (Figure S5J). 

MtAGa and MtAGb gain-of-function analyses 

Ectopic expression of results in the homeotic conversion of 
sepals and petals into carpels and stamens, respectively [55,56]. 
Other reported phenotypes include early flowering and curling of 
leaves. To investigate whether MtAGa and MtAGb proteins differ 
in their ability to induce rej)roductive organ fate we overexpressed 
these two genes in Arabidopsis. The coding sequence of these 
genes was located under the control of the CaMV 35S constitutive 
promoter. We analysed 75 independent primary transgenic (Tl) 
lines for every construct {35S::MtAGa and 35S::MtAGh plants) 
and found plants with different floral phenotypes that we classified 
as mild, medium and strong phenotype (Figure 5B-D). Transgenic 
plants with mild phenotype showed wax accumulations on the 
surface of the sepals (Figure 51) that are characteristic of the carpel 
(Figure 5H). The petals are narrow and asymmetric (Figure 5B 
and 5J) with partial homeotic conversions to stamens (Figure 5K). 
Flowers with medium phenotypes showed curved sepals (Fig- 
ure 5C) with abundant accumulation of wax (Figure 5L). Tl 
plants with medium and strong phenotypes showed near complete 
homeotic transformations of sepals to stamens (Figure 5M and 
5N). Finally, flowers with strong phenotype showed thickened and 
curved sepals (Figure 5D) that contain carpelloid structures with 
ectopic ovules and stigmatic papillae (Figure 5P). Although these 
three phenotypical classes were observed for both constructs, we 
notice a higher proportion of plants with medium/strong 
phenotypes (80.3%) in the MlAGh overexpressing lines (Fig- 
ure 5Q). In contrast, 59.5% of MlAGa overexpressing lines 
showed mild homeotic transformations. 

In summary, constitutive expression of MtAGa and MtAGb in 
Arabidopsis are able to cause equivalent homeotic alterations in 
the flowers that are essentially identical to those caused by 
constitutive i'xpr(;ssi<)n of other eu^G genes [56]. However 
MtAGb protein seems to be more efficient than MtAGa to induce 
reproductive organ fate when overexpressed in a heterologous 
system. 

Discussion 

The model legume Medicago truncatula harbors two members 
of the euAG sub-clade [MlAGa and MlAGh) and one member of 
PLENA [MtSHP) [30]. In this study, we report the functional 
characterization of the eu^G homologs from this species to 
achieve a more thorough understanding of the evolution of the 
AGAMOUS subfamily in core eudicots. 



Two euAG genes in Medicago truncatula 

The finding of two euAG duplicated genes [MtAGa and MtAGb) 
in M. truncatula, as well as in other legume species (see Figure IB) 
suggests that each pair of paralogs arose by a duplication event 
prior to the speciation of legumes. They could have been 
originated during the whole-genome duplication (WGD) event 
that predated speciation of M. truncatula and other legumes 
around 50-60 Mya [57,58]. Gene duplications allowed the 
emergence of critical changes in several MADS-box gene lineages 
that control floral organ identity. In Medicago, B-function genes 
are also duplicated: MlAP3-lihe genes were originated by an 
ancient duplication concomitant with the base of the core eudicot 
radiation [49], whereas the two duplicated P/-like genes were 
probably originated during the WGD event that also originated 
the two MtAG genes ([48] ; Roque ei at. unpublished). The changes 
of these genes after duplication may have played an important role 
in the evolution of floral morphology and ontogeny in legumes. 

The presence of two eu/lG genes is not only found in legumes. 
Previous phylogenetic studies have revealed the existence of two 
members of the eu/4G-lineage in several core eudicots [12]. This is 
the case of the paralogous CUMl and CUM2 from Cucumis 
sativus, PTAGl and PTAG2 from Populus trichocarpa, and 
GAGAl and GAGA2 from Gerhera hyhrida. These studies are 
mainly focused on the analyses of their expression patterns and 
limited functional information is available. Recendy, phylogenetic 
analyses using AG protein sequences from several Asteraceae 
species, have revealed that some of them, harbor two or three 
members within the eu^G-lineage [24]. 

MtAGa and MtAGb proteins conserx^e characteristic domains 
of the C-lineage members, as the presence of an N-terminal 
extension preceding the MADS domain [59]. The length of this 
sequence in MtAG proteins falls within the range previously 
described for other members of this clade (from 13 to 52 amino 
acids) [59] (Figure SI). However, regarding to DNA sequence, 
both genes possess six introns rather than the typical eight of the 
other y4G homologs [7,8,11,25,60,61]. This genomic organization 
is observed in both angiosperm and gymnosperm genes suggesting 
that the presence of eight introns is likely to be primitive in AG 
subfamily [12]. The genomic organization ofM. truncatula euAG 
genes could have taken place during the WGD (■\ (-nt that predates 
legume speciation followed by gene rearrangement. 

Within the highly similar gene structure ot AGAMOUS genes, it 
has been shown that the first and second introns are essential for 
their correct expression [24] . Among the identified enhancers, the 
transcription factor LFYis essential for AG activation and binds to 
CM-elements located in the second intron [52,62]. The presence of 
LFY-binding sites has been reported to be critical to establish early 
expression of C-function genes during flower development [63]. 
We identified high affinity binding sites for LFY in the first intron 
of both MtAGa and MtAGb. Although the binding site landscapes 
observed for each eu^G-paralogous gene are different, experi- 
mental analyses have demonstrated that both types of distribution 
of LFY binding sites are functionally relevant [53]. In contrast, in 
our analyses MlSHP was not found to be a likely LFY target and 
accordingly, MtSHP expression is exclusively detected in the 
ovules at later stages of floral development. Similarly, SHP or 
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Figure 5. Overexpression of MtAGa and MtAGb in Arabidopsis thaliana: floral phenotypes. Wild-type Arabidopsis flower (A). Floral 
phenotypes observed in Arabidopsis plants overexpressing MtAGa or MtAGb genes have been classified into mild (B), medium (C) and strong (D). 
Scanning electron micrographs showing the characteristic cellular types of wild-type sepals (E), petals (F), stamens (G) and carpels (H). Scanning 
electron micrographs of sepals from overexpression lines with mild phenotype (I) showing wax accumulation (arrow). Narrow petals from 
overexpressing lines (J) showing staminoid cellular types (K). Wax accumulation in sepals (L) of plants with medium and strong phenotypes. Near 
complete homeotic conversion of petal into stamen (IVI) showing characteristic cellular types from anther (N) and filament (O). Flower with strong 
phenotype showing ectopic ovules and stigmatic tissue (P) in the first whorl. Relative proportion of phenotypes observed in plants overexpressing 
MtAGa or IVltAGb genes (Q). Bars indicate: 1mm in A, B, C, and D. 
doi:1 0.1 371 /journal.pone.01 03770.g005 



STK, which play late roles in fruit and ovule development, were 
not found to be LFY targets in ChlP-seq experiments [53]. 
Together, our analyses suggest that both MtAG genes have the 



same probability to be directly activated by LFY, conferring organ 
identity on floral meristems that arise after the transition to 
reproductive development. 
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Two eoAG genes control C-function in Medicago 
truncatula 

A typical C-function knock-out phenot^pe is described as the 
conversion of stamens to petals, and a flower in place of the fourth 
whorl, which starts with a whorl of sepals [24] . 1 % of MtAGab- 
VIGS flowers showed a near complete loss-of-C-function pheno- 
type. Stamens were replaced by petals or petal-like organs and the 
central carpel was replaced by multiple petaloid or sepaloid 
concentric structur(;s. MtAGab-VIGS plants still showed some 
reproductive structures probably due to the fact that these Hnes stiU 
produce some MtAGa and MtAGb transcripts. It would be 
desirable to measure the MtAG mRNAs levels in MtAGab-VIGS 
flowers to validate this hypothesis, but sufficient plant material was 
not available. However, using stable MMG::RNAi lines we were 
able to correlate the amounts of MtAG gene products with partial 
loss-of-C-function phenotypes. Only simultaneous down-rc^gula- 
tion of both genes (below 40% of the wild-type level) produces 
mild aberrations in stamen and carpel development. These results 
support the quantitative model of AG activity, which suggests that 
the different y4G functions vary with respect to the amount of gene 
product required [47] . The floral determinacy requires the higher 
levels of AG, while low levels of AG are sufiicient for carpel and 
stamen specification [51]. According to this model, only a strong 
reduction of both MtAG transcripts could explain the indetermi- 
nacy of the floral meristem in MtAGab-VIGS flowers. In each 
single mtaga and mlagh mutants the amount of the remaining 
paralog is sufficient to confer almost all the C-function. However, 
the discrete mutant phenotypes of these plants suggest that a single 
MtAG gene is not sufficient to maintain the overall dosage of these 
subcomponents of protein complexes which are needed to specify 
the C-function in M. truncatula. 

The spatial and temporal expression profile of AGAMOUS 
subfamily genes is closely consistent with their conserved function 
and subfunctionalization [24]. The expression of AG lineage genes 
early in the florad meristem, and in the third and fourth floral 
whorls is probably the ancestral expression profile of these genes 
[24]. An overlapping conserved expression pattern for both eoAG 
and PL£-like genes is observed in Petunia, N. henihamiana and 
Antirrhinum, at early floral developmental stages. This expression 
profile correlates with their functional redundancy to promote the 
C-function of these genes in the cited species [4,11,20]. In A. 
thaliana, the only C-function factor is AG, because, although the 
PLE-like proteins (SHP proteins) fully maintain a similar potential 
activity, they are not physically present in the meristem and 
primordia cells at the appropriate time. Their functions are rather 
restricted to regulate specific tissue types after carpel identity has 
been established. It remains possible that SHP genes redundantiy 
regulate carpel identity with AG because they are expressed in 
early carpel primordia. The Arabidopsis AG and SHP genes are a 
clear example of subfunctionalization as a result of a change in 
expression pattern [64]. M. truncatula could be considered an 
illustrative example, where the expression pattern of the members 
of eu^G or PLENA Kneages would be very well correlated with 
their function. The two evAG genes [MtAGa and MtAGb) are 
expressed early during floral development, having a critical role in 
the specification of reproductive organs and floral determinacy. 
MtSHP expression appears late during floral development after all 
organs have been specified, suggesting it does not contribute 
significandy to the diflFerent components of the C-function (stamen 
and carpel identity, and determinacy). A role for this gene in fruit 
morphology and valve margin lignification has been suggested 
[30]. However, we cannot exclude a function for MtSHP in ovule 
development, associated with its expression in this tissue during 
floral development. 



It seems that changes in coding sequences are not the main 
cause of the functional divergence between the majority of the 
eu^G and PLENA-like members. In most cases, the ectopic 
expression of each sub-clade member gave similar results 
[4,23,55,65,66,67]. This demonstrates that both types of proteins 
have conserved their biochemical interactions. However, an 
exception to the rule has been observed, for example, in 
Antirrhinum, where the presence of an additional glutamine in 
the K3 helix of FAR protein, causes its inability to induce carpel 
identity in the first whorl [68,69]. These data suggest that the 
functional diversification oi Antirrhinum C-lineage genes largely 
relies on changes in the protein ability to interact rather than on 
changes in their expression patterns. The ectopic expression of M. 
truncatula eu^G- and PLENA-like genes {MtAGa, MtAGb and 
MtSHP; this work; [20]) in Arabidopsis was able to cause 
homeotic alterations in the flowers that are essentially identical 
to those caused by the constitutive expression of other cuA G- and 
PLENA-like genes. However, we noticed that although there are 
not qualitative differences in the floral phenotypes between 
33S::MtAGa and 35S::MtAGb transgenic plants, they differ 
quantitatively, indicating that MtAGb protein could be more 
active in Arabidopsis than MtAGa. 

Overall, our results show that MtAG genes could have the same 
probability to be regulated by LFY at early stages of floral 
development, while the prediction of binding sites for LFY in the 
first intron of MtSHP indicated a low probability to be regulated 
by LFY. Both MtAG genes showed nearly identical spatial 
expression patterns, which are similar to those described for C- 
function genes from the euAG subclade such as ^G and PLE. In 
contrast, MtSHP is exclusively detected in the ovules at late stages 
of floral development. It is not physicaUy present in the meristem 
and primordia ceUs at the appropriate time to participate in the C- 
function. Therefore, in M. truncatula, the C-function is redun- 
dantly encoded by the two eu^G-genes. Finally, the overexpres- 
sion oi MtAGa or MtAGb in Arabidopsis, gave similar results to 
those observed for the ectopic expression of other eu/4G-genes. M. 
truncatula tuAG- and PLE AM -like paralogs could have sub- 
functionalized concomitantly with the differential spatial and 
temporal expression pattern of the ancestral gene lineage. 
Furthermore, MtAG paralogs are maintained in the genome 
consistent with the gene balance hypothesis, which predicts that 
the fate of duplicated genes largely depends on the maintenance of 
the stoichiometric balance among members of the macromolec- 
ular complex [70,71,72,73]. 

AG subfamily has been shaped by a complex history of gene 
duplications leading to various degrees of redundancy and/or 
subfunctionalization of the C-function in individual gene copies. 
Despite these differences, the control of the development of 
reproductive organs is strongly conserved across eudicots. We 
provide insight on the evolution of the A G subfamily in a new core 
eudicot species, highlighting the dynamic nature of functional 
evolution following gene duplication in this MADS-box gene 
subfamily. 

Supporting Information 

Figure SI Sequence comparison of MtAGa, MtAGb and 
related MADS-box proteins. Identical amino acid residues are 
shaded in black, and similar amino acid residues in grey. The 
MADS and K domains, and I and C-terminal regions are marked 
with arrows. The AG motifs I and II are underlined with 
continuous and discontinuous lines, respectively. A dotted line 
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before the MADS-box indicates the N-terminal sequence 

characteristic to some C-hneage genes. 

(TIF) 

Figure S2 Southern blot analysis o£ MtAGa and MtAGb 
genes and localization of the MtAGb gene on the Af. 
truncatula physical map. (A) Southern blot analysis oi MtAGa 

and MtAGb. The length of the DNA markers (in Kb) is shown on 
the left margin. (B) MtAGb sequence is anchored in BAG 
AC15346(), on the chromosome 8. Pseudochromosomes are 
represented by the succession of BACs (vertical, grey thick lines) 
separated by gaps (black, horizontal lines). Not aU gaps are visible 
at this resolution. Unmapped BACs are shown on the right as 
MtChrO. MPMl/2: molecular weight markers. 
(TIF) 

Figure S3 Molecular characterization of Tntl insertion 
line in the mtagb locus. Genotyping of 20 plants from the Tntl 
insertion population of NF4908 line. (A) PGR results using AGb- 
F/Tntl-F primers. (B) PGR results using AGb-F/AGb-R genomic 
primers. Homozygous mutant plants are indicated with bold, 
numbers. MPM: molecular weight markers; WT: wild-type 
control; HgO: distilled water. 
(TIF) 

Figure S4 qRT-PCR expression analyses of loss-of- 
fiinction plants. Relative expression of MtAGa or MtAGb 

genes in flower buds of mlaga and mtagb mutants. The height of 
the bars for a given gene indicates differences in relative expression 
levels in floral buds. The wild-type expression value was set to 

I. 00, and lower values are plotted relative to this value. 
(TIF) 

Figure S5 Floral phenotype and expression analyses of 
35S::RNAi-Mfc4G plants. (A) Wild-type M. truncatula flower. 
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